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Abstract 

We consider three problems for Gabor frames that have recently 
received much attention. The first problem concerns the approxima- 
tion of dual Gabor frames in I^QR) by finite-dimensional methods. 
Utilizing Wexler-Raz type duality relations we derive a method to ap- 
proximate the dual Gabor frame, that is much simpler than previously 
proposed techniques. Furthermore it enables us to give estimates for 
the approximation rate when the dimension of the finite model ap- 
proaches infinity. The second problem concerns the relation between 
the decay of the window function g and its dual 7. Based on results 
on commutative Banach algebras and Laurent operators we derive a 
general condition under which the dual 7 inherits the decay properties 
of g. The third problem concerns the design of pulse shapes for or- 
thogonal frequency division multiplex (OFDM) systems for time- and 
frequency dispersive channels. In particular, we provide a theoretical 
foundation for a recently proposed algorithm to construct orthogonal 
transmission functions that are well localized in the time-frequency 
plane. 
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1 Introduction 



Gabor systems play an important role in signal processing and digital commu- 
nication. In filter bank theory they are known under the name oversampled 
modulated filter banks ||, in wireline communications they correspond to the 
concept of discrete multitone transmultiplexing, and in wireless communi- 
cations they are (implicitly) used in orthogonal frequency division multiple 
access systems [0, |25|, [3J . 

A Gabor system consists of functions of the form 

gna, mb (t) = e 2mmb g(t -no), n,meZ, a, b G R, (1) 

where g G L 2 (R) is - depending on the context - called window, atom, or 
pulse shape. The parameters a and b represent the time-shift and frequency- 
shift, respectively 

We say that (g,a,b) generates a Gabor frame for L2(R) for given shift 
parameters a, b if there exist constants (frame bounds) A, B > such that 

^ll/H 2 < £ \(f^na, mb }\ 2 <B\\ff, (2) 

for any / G L 2 (R). 

The analysis operator T is defined as 

T-.fE L 2 (R) -> T/ = {{f,gna,mb)}n,meZ, (3) 

and the synthesis operator, which happens to be the adjoint of T is 

T* : c G £ 2 (Z x Z) -> T*c = £ c ™^™, m 6. ( 4 ) 

The Gabor frame operator is defined by 

Sf = £ (fi9na,mb)g na,mbi f e L 2 (R), (5) 

and satisfies 

IA<S <IB 

where / is the identity operator on L 2 (R). Of course S = T*T. 
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If (g,a,b) establishes a Gabor frame for Lg(^) then any / in Lz^M.) can 
be represented as 



/ ^ {f j 1n,m) 9na,mb ^ (/; 9na,rab)^in,mi (6) 
n,m£Z n,m£Z 

where the dim/ frame {~i n ,m} is given by 

7n,m = e 2wim6 7(t - no), n,m eZ; a,b e R, 

with 7 = S^ x g. In general there are many functions generating dual frames 
that satisfy relation (0). The "canonical" dual window 7 has several nice 
properties. One of them is that it has minimal l^-norm among all dual 
functions. In this paper we concentrate on the canonical dual window and 
henceforth simply talk about the dual window and the dual frame. For more 
details about the properties of Gabor frames and their duals the reader is 



referred to [B, WM . 



The rest of the paper is organized as follows. In Section |2| we analyze 
the problem of approximating the dual window by using finite-dimensional 
methods. Based on Wexler-Raz type duality relations we derive a method 
to approximate the dual Gabor frame, that is much simpler than previously 
proposed techniques. Furthermore we show that for windows with exponen- 
tial decay in time and frequency the proposed approach yields an exponential 
approximation rate when the dimension of the finite model approaches in- 
finity. In Section [5] we dig deeper into decay properties of g and its dual 
7. Based on results on commutative Banach algebras and Laurent operators 
we derive a general condition on the decay of g which guarantees that the 
dual 7 inherits these decay properties. Finally in Section [| we demonstrate 
the relevance of the results derived in Section |^ for wireless communications. 
In particular, we provide a theoretical foundation for a recently proposed 
algorithm to construct an orthogonal frequency division multiplex (OFDM) 
system with good time-frequency localization properties. 

Before we proceed we introduce a few notations used throughout the 
paper. 

The Fourier transform of a function / is given by 



f(u) = / f(t)e-^dt 
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The short time Fourier transform (STFT) of / with respect to the (suffi- 
ciently nice) window g is 



(V g f)(t,u)= J f(x)g(x-t)e- 2 ™"dx. (7) 

— oo 

A locally integrable function w is called weight function, if w is positive 
and submultiplicative, i.e., if w(t) > and w{t\ + 1 2 ) < w{t\)w{t2). 
The space Z, ljM ,(IR) consists of all functions / with 



\f(t)\w(t)dt < oo, 



where w is a weight function. Similarly £i )W (Z) consists of all sequences 
x = {x k } ke z with 



\x k \w(k) < oo. 



k=oo 



It is convenient to define following spaces: 

C M = {/ G L 2 (R) : 3a,/5 such that |/(t)| = ,V|t| £ [a,/?]}. 

g A = |/ e L 2 (R) : 3 A > 0,c> such that \f(t)\ < ce" A|t| ,V|t|}. 
Q s = {/ e L 2 (R) : 3s > l,c> such that |/(t)| < c(l + ,V|t|}. 



Finally, the Moore-Penrose inverse |12fl of a bounded operator T is de- 
noted by T + . 

2 Gabor frames, finite sections and the dual- 
ity condition 

The theoretical concepts for Gabor analysis are usually developed for infinite- 
dimensional function spaces, notably for L^iM.) or whereas all nu- 
merical implementations have to be done within a finite-dimensional frame- 
work. The connection between Gabor systems on ^(Z) an d several finite- 



dimensional models has been clarified in [29, 31]. In this section we extend 
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these results to Gabor frames for L 2 (M.). Different techniques than those 
used for £2(^1) are required for X 2 (K). 

Our approach relies on a remarkable property of Gabor frames, whose 
discovery has its origin in a paper by Wexler and Raz [32]] , their result was 



later made precise and extended by Janssen j24|, Ron and Shen |p8fl , and 



Daubechies, H. Landau, and Z. Landau 111 



For given g,a,b we define the operator H by 

Hf-={(f,gk/b,i/a)}k,i£Z, f e L 2 

where 

g k/b ,i /a (t)=g(t-k/b)e 2 ^ a . 

The adjoint H* of H is 



H*c = 2J c k) ig k /b,i/ a , c = {c k j} k j e z E £ 2 (Z 2 ^ 

k,l 



We identify .f/TJ* with its matrix representation, given by 

HH* = {(g k '/b,l'/a, 9k/b,l/a)}k,l;k',l'EZ- 

The Wexler-Raz biorthogonality relations imply that the dual window 
satisfies O : 
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H*{HH*)- l a = abJ2 [(MT'Wao (8) 



where cr = {(ab)8 k0 Si }k,iez- 
Furthermore, there holds: 



Theorem 2.1 ( ||28|| ) (g, 1/b, 1/a) is a Riesz basis for its closed linear span 
if and only if (g, a, b) is a frame for ^(M) . 

Formula (|8]) and Theorem |2.1| are the main ingredients for our approach 
to approximate the dual window 7 using a finite-dimensional model. 

Now, for x E ^ 2 (Z 2 ) and n E N define the orthogonal projections P n by 



k.l 



x k: i if max{|fc|, |/|} < n, 
else. 
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We identify the image of P n with the (2n + l) 2 -dimensional space 
and write 

H n f := P n Hf = {(f,gk/b,l/a)}\k\,\l\<n- 

The matrix 

H n H* = (HH*) n = P n HH*P n = {(gk'/b,l'/a, 9k/b,l/a)}\k\,\l\;\k'\,\l'\<n- 

is a finite section of the infinite-dimensional Gram matrix HH*. 

We say that the finite section method is applicable to HH*, if, beginning 
with some neN, for each y G range (H) the equation 

H* n [H n H*]- l x (n) = PnV 

has a unique solution x^ 71 ' G Im(P n ) and as n — >■ oo the vectors tend to 
the solution of HH*x = y. 
We set 
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(n) 



for n = 0, 1, 2, . . . . 
There holds: 



Theorem 2.2 Lei (g,a,b) generate a Gabor frame for L2 
the dual window. Then 



and let 7 be 



17-7 



(n)| 







for n — > 00. 



We need following result for the proof of Theorem |272 . 

Lemma 2.3 (Lemma of Kantorovich, ||27|| ) Let {K n } ne ^ be a sequence 
of invertible operators on a Banach space B, and assume that the sequence 
is uniformly bounded above and below, i.e, that there exist constants < 
Ci, C2 < 00 such that 



Ci 



< \\K n f\\ B <C 2 



for all nGN and f G B. 



(9) 



T/ien {K n } converges in the strong operator topology to an invertible operator 
if and only if the same is true for {K~ l } n , and then 



{limiQ- 1 = lim{^- 1 } 
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Proof of Theorem 

Theorem ^7i] implies that (g,l/b,l/a) is a Riesz basis for its closed span 
with Riesz bounds A, B. Any finite subset of a Riesz basis is again a Riesz 
basis for its closed span, cf. Chapter 1 in |33f1 . Thus for {gk/b,i/a}\k\,\i\<n, n E N 
there exist constants A n , B n such that 

A < A n < H n H* n < B n < B. (10) 



Hence we can apply Lemma [2.3| to establish that the finite section method 
is applicable to HH*. 

Finally note that H n obviously converges to H pointwise for n — > oo, thus 
H'yW — i|| — > for n — > oo. □ 

It is interesting to ask if we can give some estimate on the rate of approx- 
imation of the finite section method in case the window g satisfies certain 
decay conditions in time and/or frequency. In the following theorem we con- 
centrate on windows with exponential decay in time and frequency domain, 
but it will become clear from the proof that similar results may be obtained 
for other types of decay. 

Theorem 2.4 Let (g,a,b) be a Gabor frame for -L 2 (R) and assume that there 
exist constants C,D > such that 

\g(t)\ < Ce-W, and \g(u)\ < De~ x ^ , (11) 

then there exists a X' < X and a constant C depending on the frame bounds 
and on X' , but independent of n such that 

||7-7 (n) || < C'e- X ' n . 

Proof: First we show that ([TTj) implies that the entries of HH* satisfy 

{HH*) kl . m , < L> 2 e- Al (l fc -*'l + |W'l) (12) 

for some Ai < A and some constant D 2 . 
It is clear that 

\(HH*)kl,k'l'\ = \ (9k/b,l/a,9k'/b,l'/a}\ = I (9, 9k-k',l-l') \ (13) 

and 

(9,9k/b,i/a) = {V g g){k/b,l/a) k,leZ. (14) 
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Now 

+00 



(V g g)(t,u)\= g(x)g(x-t)e- 2 ^ dx< \g(x)\\g(x - t)\dx 



—00 —00 
+00 +00 



< C 2 J e- A|x| e- A|x_t| da: < C 2 J e-^e'^'^dx (15) 



for some Ai < A. Set e — A — Ai, then 



-00 +00 
C 2 I e- xlxl e~ Xllx - tl dx < C 2 [ e~ A|:r| e- Al(l * H:E| W 



-00 —00 
+00 



< C 2 / e- £|x| e- Al|t| rfx « C 2 e- Al|t| / e~ £lxl d 



=>\{V g g){t,u>)\ <C ie - Al W (16) 

for some constant Ci depending on Ai. 
Similarly 

|(V s <?)(t,u;)| = |(V^)K-t)| < Ae-^H. (17) 
By combining flTBD and ( |I7|) and taking the square root we get 



\{V a g){t,u)\<y/C^Dle-^^. (18) 

This together with equations fll3D and (|14|) yields (|T2l). 
Now consider 

||7-7 {n) || < \\H*(HH*y 1 o — H*{HH*y l H n H* n (H n Hn)~ l (? i:n ^\\ + 
\\H*(HH*Y 1 H n H^ l (H n H^y l a^ - H*(HH*)- l HH* n (H n H* n )- l ^ n) \\ 
< \\H*(HH*)~ 1 ((T - (tV>)\\ +A- 1 \\(H n H* n - HH^H^y^W 







= A- l \\{H n H: - HH* n ){H n Hl)- x a^\\. 

Note that (H„H* — HH*) is a matrix that has a finite number of columns 
and a biinfinite number of rows, of which the rows with index \kl\ < n are 
zero. 
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By (|T2|) the entries of H n H^ and the nonzero rows of (H n H* - HH^) 
decay exponentially off the diagonal. Relation fllCiD implies 

cond(H n H*) < cond(HH*) for n = 0, 1, . . . . (19) 



Using (19) and Proposition 2 in [23] it follows that there exists a A2 and a 
constant C 2 depending on A2 and on the condition number of HH*, such 
that 

\[(H n H:)-% l;0>0 \ <C 2 e~ x ^ + W, 

and the same is true for v (n) := (H n H* - H H*)[(H n H*)-\i }0fi . Now it is 
easy to see that there exist A' and C such that 

|| V (»)||= \ v m\ 2 < C'e~ x ' n . 

\k\,\l\>n 

□ 

For more results on the connection between the decay of a function and 
its short time Fourier transform the reader is referred to [^(J . 

To fully appreciate the simplicity and advantages of the proposed ap- 
proach we consider for comparison the following method for determining the 
"dual" expansion coefficients (f,j n ,m)- 

The function / 6 L 2 (M.) can be expressed as 

/ = ^ Cn,m9na,mb (^0) 

k,iez 

where (TT*) + c = d, (21) 

with d = {d n , m } = {(/, g n a,mb)} = Tf. We can identify TT* with its Gram 
matrix representation (TT*) m>n = (g n > a ,m'b, g na ,mb) for m,n,m',n f e Z. 

Setting T n T* = P n TT*P n and d^ = P n d, we obtain the n-th approxima- 
tion to c by solving 

T n T n *c (n) = d {n) . (22) 

Unfortunately the (generalized) inverse of T n T* is not bounded for n — > oo, 
although \\(TT*)+\\ < A" 1 . In fact, \\(T n T*)+\\^ oo for n -> oo. Hence c< n > 
does not converge to c and without further modifications this approach does 
not lead to an approximation of /. 
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Instead of computing (T n T*) + we can compute a regularized inverse via a 
truncated singular value decomposition by setting the singular values of TT* 
below a certain threshold r to zero. Let (T n T*) T ' + denote this regularized 



inverse. It is shown in [^2 that (T n T*) T ' + converges strongly to (TT*) + if we 
allow the threshold parameter to vary for each n. 

In order to use this approach for practical purposes we need good es- 
timates for the optimal threshold r n . Assuming a numerical precision of 5 
of the data and setting B n = \\T n T*\\, it is shown in that r n can be 
estimated by 

'S\ £+i 



< B 



n 



Here - without going into details about regularization theory - p can be seen 
as "smoothness parameter" [T2j , the standard setting for p in regularization 



theory is p = 1 or p = 2. Thus for large n we get 

f G [B(5/2)*,B6*] 

where B is the upper frame bound. Good estimates for the upper frame 
bound are important to apply this method in practice. 

All this extra effort is not necessary when using Theorem |2.2| . A different 
approach has been proposed in ||, which however requires substantially more 
effort compared to the Wexler-Raz based approach in this section. 



Laurent operators and decay of dual Gabor 
frames 



A natural question for Gabor frames - also in spite of Theorem |3.4| is the 
following. Given a window g with certain decay properties in time and/or 
frequency, does its dual 7 have the same decay properties? This question 
is not only of interest from a theoretical viewpoint, but has a number of 
practical implications, e.g., see Section |. 

We briefly summarize a few important results. We assume in the sequel 
that {g na ,mb} constitutes a Gabor frame for 2v 2 (lR) or £2^) ■ 
(i) If g G C[a,/3], then 7 G C^^] only in very special cases 0. In general 7 



is no longer compactly supported, but has exponential decay (see |29f for a 
proof in £2(^)5 the result can be easily extended to 2v 2 (IR)). Hence in this 
case g and 7 do not belong to the same type of space. 
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(ii) If g G £\ then 7 G £\ 1 , but in general with a different exponent ji < 
7 ]29], ||. Thus g and 7 have the same type of decay, but do not belong to 
the same space, we lose some quality of decay.Q 

(hi) If g G Q s , then 7 G Q s (see |29| for a proof for ^2(2^), the result can be 
easily generalized to -L 2 (M) using the same approach as in |§). In this case 
g and 7 actually belong to the same space Q s . 

From an algebraic point of view case (hi) is the most appealing one. Can 
we find a (simple) condition on the decay of g that implies that the dual 7 
belongs to the same space, similar to case (hi)? In this section we will give 
an exhaustive answer to this question. 

We need some preparation before we proceed. 

Definition 3.1 Letw(t) be a continuous weight function onf . The (weighted) 
Wiener Algebra A w is the Banach space of absolutely convergent Fourier se- 
ries of period 1 (cf. f$Bj), i.e., f G A w if 



f(t) = Yl Cke 



2wikt 



(23) 



with 



The norm on A w is 



\ck\w{k) < 00. 



(24) 



k=— 00 



00 

\ck\w k . 

k=—oo 



It follows from Chapter 19.4 of |17 that A w is a Banach algebra under 
pointwise multiplication. A w can be identified with the space of all sequences 
c = {cfc} which are in £i iW . 

Definition |3.1| can be extended to matrix-valued functions in a straight- 



forward manner. Put 



2nit\ 



'fn(e 2mt ) 
f m i{e 2mt ) 



fim{e 



2mt\ 



fn 



^2mt\ 



(25) 



1 Note that there do exist Gabor frames with Gaussian decay, that have non-canonical 
duals with Gaussian decay, see Example 3.10 in S. 
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and set 



i 

A k = J $(e 
o 



27rit ) e - 2nit dt. 



(26) 



Let the weight function w(t) act on 
Wiener algebra A™ if 



Then $ belongs to the matrix 



|v4fc|u;(A;) < oo. 



(27) 



k= — oo 



$ is unitarily equivalent to the block Laurent operator L whose matrix 
representation (with respect to the standard basis) is given by 



A 




A- 2 


A 1 


A 




A 2 


A 1 


^0 



Here Aq denotes the (0, 0) entry which acts on the 0-th coordinate space. 
$ is also called the defining function of the block Laurent operator L. By a 
slight abuse of notation we also write L = [Aw]fc/ = _ tX) where Am = Ak-i- 
We define W™ as the space consisting of all block Laurent operators whose 



m x m blocks A^ satisfy (|27|). If m = 1 we simply write W w and L reduces 
to a scalar-valued Laurent operator. 

Theorem 3.2 Let L £ W w fre self-adjoint and positive definite. Assume that 
the weight function w satisfies 

1 



lim - 



1 and lim n y/w{n) = 1, 



(28) 



then L' 1 G W„ 



Proof: Set L = |/l;w]u=-oo and L 1 = [-B w ]^ = _ oc . Since L is positive 
definite we have 



(29) 



k=— oo 
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and by the properties of Laurent operators fT8" 



oo 

2niku 



k=—oo 



The property L E W w is equivalent to f E A w . By Theorem 2 on page 24 
in |17| an element of A w has an inverse in A w if it is not contained in a 
maximal ideal. Any maximal ideal of A w consists of elements of the form 
(cf. Chapter 19.4 in Q) 



where £ = pe with 



and 



Pi< P< P2 



pi = lim — — = and p2 = lim n y/w^. 

Due to assumption (p^) we get p\ = p 2 = 1, hence p = 1. Thus a necessary 
and sufficient condition for an element in A w to be not contained in a maximal 
ideal is Ylk dke 2mkuJ 7^ for all uj. By assumption L is positive definite, hence 
f(u) = J2k A k e 27rikuj > for all u, consequently L' 1 E W w . □ 
The proof of Theorem |3.2| is essentially based on results of Gelfand, 



Raikov, and Shilov [17|]. A crucial role plays condition (|28D , to which we 
will henceforth refer as the Gelfand- Raikov- Shilov condition (GRS-condition 
for short) . 

Now let L be a block Laurent operator with defining function $ = 
[fij]™j=i £ A™. If L is hermitian, positive definite then 

det([^(A)]- =1 )>0, |A| = 1. 

Hence we can apply Theorem 8.1 on page 830 in [113] to extend Theorem W?2 



to block Laurent operators and obtain the following 
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Corollary 3.3 Let L G be self-adjoint and positive definite. Assume 
that the weight function w defined on M m satisfies the GRS- condition, i.e., 

lim - 1 = 1 and lim n y/w(n) = 1, (30) 



t/ien L- 1 G W™. 

It is clear that in the derivations above we can replace 1-periodic functions 
by a-periodic functions (a G Q). Similarly, sequences with indices in Z can 
be replaced by sequences indexed by Z Q (i.e. by indices of the form ka). 

Now we are ready to prove a general result on the decay properties of 
Gabor windows g and their duals 7. 

Theorem 3.4 Let (g, a, b) generate a Gabor frame for -L 2 (M) with ah = - G 
Q wt/i relative prime integers p and q. Let g G .Li^QR) where w satisfies the 
GRS- condition. Then 7 G -L lu ,(IR). 



Proof: The assumption g G Li tW (M) implies J^^l.^ \g(k)\w(k) < 00 and 
£^_oc \g{k/b)\w{k/b) < 00 for 6e N. Denote 



1 k I 

Gkl ^ = lJ2 9(t~ra- -)g*(t -ra--), k,leZ. (31) 

r=— 00 

G(t) is periodic in t with period a. Since afe = | we also get Gkiit) = 
Gk+p,i+ p (t), cf. Hl3| ) Theorem 3.2. In words, G(t) is a block Laurent operator. 
The frame property implies that G(t) is hermitian positive definite and 

AI < G(t) < BI. 

The submultiplicativity of the function w implies that the spaces JDi )W (R), 
£i >w (Z), and £i jW (Zi) are Banach algebras under convolution (cf. |I7[ ). Hence 
{G ofe (t)} feeZ G £i, w (Z?), and G(t) G W™. By Corollary [T3| 



G _1 (t) G W™. (32) 

The dual 7 can be expressed as (see e.g. ||) 

00 , 
7(f) = E [G-\t)] ok g(t--). (33) 

fc = — OO 

This equation together with ( |32"D and the assumption (7 G Li )M (M) yields that 
7 G JWR). □ 
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Remark: (i) The idea to exploit the Laurent operator property of the Gabor 
frame operator in the context of window decay is not new. It has been used 
by Feichtinger and Grochenig |H| in connection with modulation spaces and 
polynomial weights. 

(ii) It is easy to reformulate Theorem |3.4j for windows whose decay properties 
are given in the frequency domain. We leave this modification to the reader. 

Observe that condition (p8|) is satisfied e.g. for w(x) = (1 + \x\) s and 
w(x) = e A ' :c ' 7 ,7 < 1, but not for w(x) = e A ' a '. This is why we have to use 
a smaller exponent Ai < A for the dual to bound the decay in case of an 
exponentially decaying window. 

A detailed discussion on the connection between decay properties of win- 
dows and modulation spaces can be found in [ITU . 



4 Orthogonal frequency division multiplex- 
ing and Gabor systems 

Orthogonal frequency division multiplexing (OFDM) has attracted a great 



deal of attention as an efficient technology for wireless data transmission |34 
Among others it is currently used in the European digital audio broadcasting 
standard ||. In a wireline environment OFDM is known under the name 
discrete multitone transmultiplexing (DMT). 

The basic idea of OFDM is to divide the available spectrum into several 
subchannels (subcarriers). A baseband OFDM system is schematically rep- 
resented in Figure [l|. 

Transmitter: Assuming N subcarriers, a bandwidth of W Hz, symbol length 
of T seconds, and carrier separation F := W/N, the transmitter of a general 
OFDM system uses the following waveforms 

^(t) = tfj(t)e 2nitlF , keZ,l = 0,...,N-l. (34) 

The transmitted baseband signal for OFDM symbol number k is 

N-l 

s k(t) = ^ c M ipi(t - kT), (35) 

1=0 

where c^q, c^i, . . . , jy-i are the complex- valued information bearing coeffi- 
cients (data symbols). Assuming an infinite sequence of OFDM symbols is 
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Transmitter 



).k ¥ (0 



u 



3 N-i,i k ¥ N _,(t) 




Channel 




(Hs)(t) 


H 









Receiver (k+l)T 



Figure 1: A general baseband OFDM system configuration 



transmitted, the output from the transmitter is a superposition of individual 
OFDM symbols 



oo N-l 



(36) 



k= — oo 



k= — oo (=0 



Receiver: The OFDM receiver consists of a matched filter bank {&} of 
similar structure as the transmitter waveforms, i.e., 



<f) kl (t) = (j)(t - kT)e lmUt , k e Z, I = 0, . . . , TV - 1. 



(37) 



The transmitted data are recovered by projecting the received signal r = 
Hs + v (where v represents additive white Gaussian noise, AWGN for short) 
onto the functions fai, i.e., 



Cki = (r,(j>ki)- 



(38) 



In the standard OFDM setup the functions ipki are designed to be mutually 
orthogonal. In this case <f) = ip. The situation where the sets {ipki} and {4>ki} 
are biorthogonal is referred to as biorthogonal frequency division multiplexing 
(BFDM). 
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It is useful to have the following two OFDM setups at hand. 
Continuous-time model: Since s(t) is a continuous signal, it is convenient 
to consider a continuous-time transmission set that is infinite in both indices 
k and I, which means we set Cm — if I < or I > N— 1. Using the viewpoint 
it is clear that OFDM can be interpreted as an orthonormal Gabor system 
in X 2 (R). For a theoretical analysis it is often advantageous to work with 
this time-continuous model. 

Discrete-time model: In practice the OFDM system is digitally imple- 
mented. At the transmitter the signal s is passed through a digital-to-analog 
converter to obtain a continuous signal. At the receiver an analog-to-digital 
converter transforms the received signal back into a discrete-time signal. 
Therefore it is useful to consider the following discrete-time setup for OFDM. 

Let ip G £ 2 (Z) and set tfj k ,i( n ) = ~ kT)e 2nilF for k,n G Z, I = 
0, . . . , N - 1 and F = W/N. Then the OFDM signal becomes 

oo N-l 

s{n) = ^2 c ^Ki- ( 39 ) 

k=— oo 1=0 

It is clear that {ipk,i} coincides with a discrete-time Gabor system. 

4.1 Wireless channels and time-frequency localization 

In the ideal case the channel H does not introduce any distortion, and the 
data Cki can be recovered exactly. However in practice wireless channels 
introduce time dispersion as well as frequency dispersion (in addition to the 
usual channel noise). The time dispersion is caused by multipath propagation 
and can lead to intersymbol interference (ISI). Frequency dispersion of the 
mobile radio channel is due to the Doppler effect and can cause interchannel 
interference (ICI). The distortion resulting from channel dispersion depends 
crucially on the time-frequency localization of the transmitter pulse shapes 
{4>m}- Robustness against doubly dispersive channels can be achieved by 
pulse shapes with good time-frequency localization. 

An optimum OFDM system in case of doubly dispersive channels would 
consist of orthogonal basis functions with TF = 1, such that the ipk,i are well 
localized in time and frequency. The condition TF = 1 (critical sampling) en- 
sures maximal spectral efficiency of the transmission system. Unfortunately, 
such a system cannot exist due to the Balian-Low theorem ||. 



17 



Therefore other approaches (cyclic prefix, pulse-shaping, BFDM) have 
been proposed for doubly dispersive channels |34|, ^Ij, ^5], Q. These ap- 
proaches can be interpreted in the time-frequency plane as using an un- 
dersampled grid, i.e., TF > 1, which results in a set of basis functions that 
is incomplete in ^(K.) (or £2^), respectively). Although the choice TF > 1 
leads to a loss in the capacity of the transmission system (determined by 
the undersampling rate), it is usually an acceptable price to pay to mitigate 
interference. 

We know from Gabor theory that there exist Gabor systems with good 
time-frequency localization for TF > 1 (take for instance TF = 2 and a 
Gaussian window). In the present context of OFDM we are particularly 
interested in the construction of orthonormal incomplete Gabor systems with 
good time-frequency localization properties. 

The following result ties tight Gabor frames for Z^O^) to orthonormal Ga- 



bor bases for subspaces of L 2 (M.). It is a simple consequence of Theorem |27T 



Corollary 4.1 ( ||28|| ) Assume \\g\\ = 1. (g,l/b,l/a) is an orthonormal ba- 
sis for its closed linear span if and only if (g, a, b) is a tight frame for ■ 

Now, let {g n a,mb} be a frame for £-2 (JR.) with frame operator S, then a 
standard way to construct a tight frame is the following. Compute 

9 = S-*g (40) 

then the set {g na ,mb} is a tight frame for X 2 (J^) and by Corollary f|4] (after 
normalizing g) {gk/b,i/a} is an orthonormal basis for its closed linear span 
(equivalently we can apply this orthogonalization procedure directly to the 
set {gk/b,i/a})- Our OFDM system is now given by setting ip := g, T = 1/6, 
F = 1/a, or equivalently, 

Ipkl ■= 9k/b,l/a- (41) 

In order to obtain an OFDM system that is well localized in the time- 
frequency plane, it seems natural to start with a function g with good time- 
frequency localization such that (g,a,b) generates a frame for L%(M.), apply 
the orthogonalization procedure ( fiff ) and hope that g = S~^g and whence ip 
inherit these localization properties. This approach is not new, it has already 
been considered in [[25|, Q], however with very different conclusions. 
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It is stated in [^5] that "such an orthogonalization of the pulses, however, 
is not desirable . . . because the good time-frequency localization of the pulses 
is destroyed by such a transformation" . In contrast in |1[ it is claimed that 
the orthogonalization procedure fl40|) u in practice . . . starting from a well- 
localized initial filter . . . yields well-localized orthogonal filters" . 

Who is right and who is wrong? The answer is: 

Both and nobody - it depends! 
The first step to this answer is the following result. 

Corollary 4.2 (i) Let g,g G L\ tW and let (g,a,b) generate a Gabor frame for 
L 2 (R). Set g = S~^g. If the weight function w satisfies the GRS- condition, 
then g,g G L^ w . 

(ii) Let g,g G l\ yW and let (g,a,b) generate a Gabor frame for ^(Z) . Set 7 = 
S~ l g and g = S~^g. If the weight function w satisfies the GRS -condition, 
then 7, 7 G £\ )W and g,g G t\^ w . 



Proof: (i) The proof of this result is similar to the proof of Theorem pO . 



We only indicate the necessary modifications. Let $ be the matrix-valued 
defining function of a hermitian positive definite Laurent operator L and de- 
note L~2 = [Cklfi^oe- It follows from the basic properties of block Laurent 
operators that > for all u and 

1 

C k = / ; 1 e~ 2nikuj duj, keZ. 

J 



The fact that L is self-adjoint positive-definite implies that L^ is in the same 
algebra as L, see p|. Thus if L G (and consequently L2 g W™) for 



weights satisfying (^), then G W™. The rest follows now by repeating 
the steps in Theorem |3.4| and using the remark following Theorem |3.4| . 
(ii) In this case the frame operator S is a block Laurent operator (see [p9fl ). 
The assumption g G t\ tW implies that S G W™. Along the same lines as above 
and in Theorem |3.4| it follows that S~ l and G and consequently 



7,7 G £i tW , similar for g,g. □ 

Thus if the window (i.e. pulse shape) g satisfies the properties of Corollary |4.2| 
the orthogonalization procedure (|40D yields a function that is in the same 
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space (in the same algebra) as the function g, and the quality of decay does 
not change. 

For a window g with exponential decay (in time and frequency) it has 



been shown that g also has exponential decay (see |2{| for a proof for £2(7*), 
this result has later been extended to L 2 (M.), cf. [|[). However, similar to the 
dual window, in general the exponent for the decay of g will be smaller than 
for g. 

Since the Gaussian is optimally localized in the time-frequency plane, in 
the sense that it minimizes the uncertainty principle, it is interesting to note 
that for g(t) = e~ nt the function g has only exponential decay in time and 
frequency, cf. 0. Thus we certainly loose some time-frequency localization 
in this case, although the resulting function g still has exponential decay. 

Thus it seems there exists a large class of windows whose time-frequency 
localization properties are not affected by applying (fTOp. This is however 
only half of the truth, since in the considerations above we have ignored any 
constants that come into play. This is acceptable from an asymptotic-analysis 
viewpoint, but not for applications. 

Consider for example the following situation. Take a window g with ex- 
ponential decay in time and frequency and assume that (g,T,F) constitutes 
a frame for TF > 1. Then for any e > there exists an N such that for 
T > N, F > N 

9 9 



< £ 

\9\\ \\9\"' 

since (TF)'y — > g for T — > 00, F — > 00 (see [Id"] for a mathematically precise 



formulation). But of course TF ^> 1 leads to an unacceptable large loss of 
capacity for OFDM, already TF > 2 seems to be prohibitive in this context. 

On the other hand, if we let TF — > 1, then g and 7 will become in- 
creasingly "ill-localized" in the time-frequency plane (since in the limit case 
TF = 1 we are confronted with the Balian-Low theorem), although from 
theory we know that g has exponential decay as long as TF > 1. However, 
as pointed out earlier, this decay involves a constant C that depends on the 
ratio of the frame bounds. Since A — > for TF —>■ 1 it follows that C — > 00. 
Thus for applications such as OFDM the exponential-decay property of 7 
and g quickly becomes meaningless for TF close to 1. If we choose a well- 



2 Due to the connection between Gabor frames and modulated filter banks 0, Corol- 
lary 4.2 also yields a characterization of the decay properties of paraunitary modulated 
filter banks. 
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localized window, there is a trade-off between increasing TF and increasing 
the frame bound ratio B/A. 

Hence a correct formulation from a practical viewpoint of the two con- 
tradicting statements above is: 

If g is well localized in time and frequency and if the frame bounds 
satisfy B/A « 1, then g is also well localized in time and fre- 
quency. 

Numerical experiments indicate that for TF = 1.3 or TF = 1.4 it is 
possible to construct OFDM basis functions with good time-frequency local- 
ization, e.g. choose the Gaussian as initial window g. A definite answer if the 
time-frequency properties of the resulting OFDM basis functions obtained in 
that way are sufficient for practical purposes is difficult, since it depends on 
the actual ISI and ICI, as well as the AWGN behavior of the channel. 

Remark: It is well-known that Wilson bases can be constructed from twofold 



oversampled tight Gabor frames [10]. In light of this fact Corollary O pro- 
vides an extension of existing results on the connection between the decay 
behavior of Gabor frames and Wilson bases. 
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